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ABSTRACT 

We investigate the expected correlation between the weak gravitational shear of distant 
galaxies and the orientation of foreground galaxies, through the use of numerical simulations. 
This shear-ellipticity correlation can mimic a cosmological weak lensing signal, and is poten- 
tially the limiting physical systematic effect for cosmology with future high-precision weak 
lensing surveys. We find that, if uncorrected, the shear-ellipticity correlation could contribute 
up to 10% of the weak lensing signal on scales up to 20 arcminutes, for lensing surveys with a 
median depth z„i = 1. The most massive foreground galaxies are expected to cause the largest 
correlations, a result also seen in the Sloan Digital Sky Survey. We find that the redshift de- 
pendence of the effect is proportional to the lensing efficiency of the foreground, and this 
offers prospects for removal to high precision, although with some model dependence. The 
contamination is characterised by a weakly negative B-mode, which can be used as a diag- 
nostic of systematic errors. We also provide more accurate predictions for a second potential 
source of error, the intrinsic alignment of nearby galaxies. This source of contamination is 
less important, however, as it can be easily removed with distance information. 

Key words: cosmology: observations - gravitational lensing - large-scale structure. 



1 INTRODUCTION 

Weak gravitational lensing by large-scale structure provides a 
test of cosmology that is based on the well understood physics 
of gravity. The lensing strength is directly related to the large- 
scale distribution of matter and thus enables constraints to be 
placed on cosmological parameters such as the matter density 
pa rameter Q.^ and the ampli t ude of the mat t er power spectrum 
(Tg^Van Waerbeke et al.' '200 1'; 'Hoekstra et al.' '2002'; 'Bacon et alj 
2003; Jarvisetal. 2003; Brown et al. 2003; Hamana et al. 2003; 
Massev et al. 2005; Rhodes et al. 2004; Van Waerbeke et al. 2005; 
Hevmans et al-2005. ; .Semboloni et al_2006; .Hoekstra et al-200Qn . 



Combined with redshift information, measurements of weak lens- 
ing are also able to constrain the equation of state of the dark energy 
'w{z) as lensing can probe the effect dark energy has on both the 
evolution of the matter power spectrum and the redshift-distance 
relationship. With the next generation of deep-field and wide-field 



* heymans ©physics. ubc.ca 

t Scottish Universities Physics Alliance 



multi-colour surveys, weak lensing has the exciting prospect of pro- 
viding constraints to distinguish between a cosmological constant 
(w = —1) and quintessence models {w{z) ^ —1). 

Lensing by large-scale structure weakly distorts images of 
background galaxies such that the observed galaxy ellipticity 



Gobs 



Cs + 7 



(1) 



where Cs is the intrinsic galaxy ellipticity, or source ellipticity, and 
7 is the weak lensing shear distortion. Lensing by large-scale struc- 
ture is coherent and hence it induces weak correlations in the ob- 
served ellipticities of galaxies such that (eoi,seobs)e > 0, where the 
average is taken over all galaxy pairs (a, b) separated by angle 9. 
From equation^ 



/a b \ / a b\ , / a b\ , / b a\ , / a b\ 

(eobsCobs> = (esCs) + (7 ej + {y e^} + {j 7 ) 



(2) 



In order to extract the weak cosmological signal {7" 7'') from 
the observed galaxy ellipticity correlation (eobsCobs)' ^H intrin- 
sic terms including eg are typically assumed to be negligible. 
There are however potential physical limitations to this assump- 
tion related to correlations between galaxy shape and the den- 
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sity field. This could arise, for example, through tidal torques 
operating during galaxy formation lHo vldll949l : |Peebleall969l : 
iHeavens & Peacock 1988). Such a correlation would introduce an 
intrinsic alignment between nearby galaxies mimicking the weak 
lensing effect, i.e if galaxies a and 6 were physic a lly clo se (e°e^) ^ 
0. In vestigated observationally fBrown et al.' ^2002'); 'Lee et alj 
J2002): Hevmans et al. (2004), hereafter H04, Mandelbaum et al. 
(2006)), analytically ( Catelan et al. 2001; Crittenden et al. 2001; 
[Lee & Pen 2001; Porciani et al. 2002; Mackev et al. 2002) and 
throu gh numerical simulations ( Heave ns et alj J200CI) . hereafter 
HRH, ICroft & Metzlen >2000i) ; iJingi t2003) . H04). it is currently 
thought that this intrinsic alignment effect is weak, at the level of 
a few percent contamination to the lensing signal from a deep sur- 
vey of median redshift Zm ~ 1. Divi sion of a redshift survey for 
lensing tomography (see for example lHiill999l) . however, signifi- 
cantly increases the contamination in the lower redshift bins (H04). 
As sho wn by [Hevmans & Heavens (2003) and King & Schneider 
yOOZ) this contamination can be effectively removed at high ac- 
curacy by optimally downweighting pairs of nearby galaxies in the 
lensing analysis (see H04 for an implementation of this method). 
Hence this physical systematic need not limit the accuracy of fu- 
ture multi-colour weak lensing surveys. 

Correlations between a galaxy's shape and its surrounding 
density field may al so introduce another more subtle effect, first 
discussed bv lHirata & Seliak (2004). Consider the image of a dis- 
tant galaxy lensed by the dark matter density surrounding a fore- 
ground galaxy. If the shape of the foreground galaxy is correlated 
with the tidal gravitational field, this could lead to a non-zero cor- 
relation^ between intrinsic foreground galaxy ellipticity es and the 
lensing shear distortion experienced by the distant galaxy 7, i.e 
if galaxy a is always more distant than galaxy h, {'y'^e^) 7^ 0. 
iHirata & Seliak ( 2004) explore two analytical models to estimate 
the magnitude of this effect finding that when galaxy ellipticity is 
quadratically related to the tidal or density field, as suggested by 
tidal torque theory, no shear-ellipticity correlation is found. When 
galaxy ellipticity is linearly related to the tidal or density field, how- 
ever, the shear-ellipticity correlation can be similar to the magni- 
tude of the weak lensing shear correlation. In this model, galaxies 
are assumed initially to be spherical and acquire ellipticity by re- 
sponding to the linear tidal field. For a spherically symmetric over- 
density, the tidal field compresses the spherical galaxy along the di- 
rection to the over-density centre. Around a spherically symmetric 
proto-void, the linear theory would stretch the galaxy along the di- 
rection to the void centre. There is some observational evidence to 
support the linear model from Truiillo et al. (2006), who find that 
the rotation axes of galaxies located near voids are preferentially 
aligned with the void surface. 

If a shear-ellipticity correlation does exist then its strength 
is expected to be dependent on the lensing efficiency of the lens- 
source pair (Hirata & Seliak 2004 ). King ( 2005) has developed 
this idea to show that using correlation function tomography and 
knowledge of the redshift dependence of each source of contamina- 
tion, one could separate the true cosmological lensing signal from 
both the intrinsic galaxy alignment contamination and the shear- 
ellipticity correlation contamination. 

The shear-ellipticity correlation effect has recently been inves- 



^ We he reafter refer to this effe ct a s a shear-ellip ticity conelation. In the 
literature i Hirata & Selia3 12.004") and'Kina' '2005^ refer to this effect as in- 
trinsic alignment-lensing interference and, Ma ndelbaum et ai., C.2QQ6.) use the 
term GI contamination. 



tigated observationally bv lMandelbaum et alj l200d) who measure 
the correlation between the intrinsic ellipticities of nearby galaxies 
and the density field that surrounds them in the Sloan Digital Sky 
Survey (SDSS). As the density field is responsible for the gravi- 
tational lensing of background galaxies, Mandelbaum et al. ( 2OO3) 
are able to estimate the shear-ellipticity correlation for any lensing 
survey. They find that for most galaxy types, the shear-ellipticity 
correlation is consistent with no signal. The shapes of the most lu- 
minous red galaxies are however found to be be strongly correlated 
with the density field. For a lensing survey with a median source 
redshift z™ ~ 1 the most luminous galaxies would then introduce 
a shear-ellipticity correlation signal that is found to be at a similar 
magnitude to the pure lensing correlation signal f rom a ACDM cos- 
mology. Averaging over the full galaxy sample. iMandelbaum et alJ 
(2006) suggest that there would be a ~ 5% contamination to the 
pure lensing signal. These results are however subject to noise (the 
results for the full galaxy sample are also consistent with zero con- 
tamination), and are based on measurements taken at very low red- 
shift and projected to the higher redshifts that are probed by lens- 
ing surveys. In this paper we therefore aim to address this issue 
through the analysis of extensive numerical dark matter simula- 
tions (Vale & White 2003), following the philosophy of HRH and 
H04 whereby galaxy ellipticity is defined from the properties of the 
parent dark matter halos. 

This paper is organised as follows. In section |2| we describe 
the N-body simulations that we use in this analysis and the different 
models used to populate dark matter halos with luminous galaxies. 
We update the results of HRH and H04 in section |3] by present- 
ing higher signal-to-noise measurements of the E and B-mode in- 
trinsic alignment correlation {e^e'^)g from our N-body simulations 
and we compare these results to the recent observational measure- 
ments of Mandelbaum et al. ( 2006) . We also investigate the effect 
of nearby voids on galaxy sh ape and compare our results to the ob- 
servational measurements of iTruiillo et alj i2006ft . In section|4]we 
present measurements of the E and B-mode shear-ellipticity corre- 
lation (7"es) and measure the redshift dependence of this signal. 
This allows us to provide fitting formulae for the shear-ellipticity 
correlation measured from different galaxy models. We discuss our 
findings in section|5]and conclude in section|6| 



2 N-BODY SIMULATIONS 

The basis of our calculation is a large N-body simulation of a 
ACDM cosmology. The simulation used 512^ particles in a pe- 
riodic cubic box 300/i^^Mpc on a side. This represents a large 
enough cosmological volume to ensure a fair sample of the Uni- 
verse, while maintaining enough mass resolution to identify galac- 
tic mass halos. The cosmological model is chosen to provide a 
reasonable fit to a wide range of observations with fim — 0.3, 
Qa = 0.7, Ho = 100/ikms~^Mpc"^ with h = 0.7, t^s/i^ = 
0.02 and a scale-invariant spectrum with as = 0.8. The simu- 
lation was started at z = 60 and evolved to z = using the 
TreePM code of lWhitd i2002l) . The full phase space distribution 
was recorded every 100/i~'^Mpc from z ~ 2 to 2; = 0. The grav- 
itational force softening was of a spline form, with a "Plummer- 
equivalent" softening length of 20 h~^kpc comoving. The particle 
mass is 1.7 x lO^^'h^^MQ allowing us to find bound halos with 
masses several times 10^^/i~^Mq. 

For each output we determine the halo catalogue by running 
a "friends-of-friends" (FoF) group finder (e.g. Davis et al. ( 1985)) 
with a linking length b — 0.15 in units of the mean inter-particle 
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spacing. This procedure partitions the particles into equivalence 
classes, by linking together all particle pairs separated by less than 
a distance b. This means that FoF halos are bounded by a surface of 
density roughly 140 times the background density. For each halo 
we computed the angular momentum and the moment of inertia 
tensor, using the centre of mass of the particles in the FoF group as 
the centre of the halo. 

A past light cone was constructed by stacking the boxes back 
along the line of sight, with boxes shifted by random amounts in 
the X- and y-directions and viewed down randomly chosen axes to 
avoid artificial periodicity. The full projected mass distribution was 
used to create convergence maps of k at a source redshift plane of 
Za ~ 1 and Zs ~ 0.5, where k is defined to be the dimensionless 
surface mass density. The gravitational shear 7 — (71,72) that a 
galaxy at Zs ~ 1, or Zs ~ 0.5 would ex perience is then created 
from the k maps rKaise r & SQuireslll993lh . The halo information 
was transformed into the same coordinate system to create a light 
cone halo distribution. 12 different realizations of the 2048^ pixel 
maps were made in all, each with a 25 square degree field-of-view. 

2.1 Galaxy Models 

From a catalogue of dark matter halos, there are two steps required 
to simulate luminous galaxies with ellipticities. Firstly we need to 
populate the halos with galaxies, for whic h, in this paper, we use the 
conditional luminosity function (CLF) of ICoorav & Milosavlievid 
yOO^, and secondly we need to assign ellipticities. We follow 
HRH and H04 in relating the ellipticities e of the galaxies to the 
ellipticity of the halo (the 'elliptical' model) or via the angular mo- 
mentum of the halo (the 'spiral' model). In this section we first 
describe the CLF halo occupation model, and then discuss the el- 
liptical, spiral and galaxy mix models that we use in the analysis 
that follows. 

The CLF, denoted by $(L|Mh, z) gives the average number 
of galaxies with luminosities between L and L + dL that reside 
in halos of mass M^ at a redshift of z. The CLF model is sepa- 
rated into terms associated with central and satellite galaxies with 
a further division into early-type and lat e-type galaxies . The free 
model parameters are then constrained in Coorav (2006) by obser- 
vations of the galaxy luminosity function and galaxy correlation 
function over a w ide range of redshi ft samples (for example the 
SDSS at 2 < 0.1 feeha vi et all2004l) and the COMBO- 17 survey 
at 0.4 < z < 0.8 lWol fetalJl2'003h ). The results of this analy- 
sis show a very convincing match between CLF model predictions 
and observational results at all redshifts and thus provide an ade- 
quate model with which to populate simulated dark matter halos 
with luminous galaxies. The average number of galaxies that re- 
side in a halo of mass Af^ at redshift z is given by N{Mh, z) = 
J dL^[L\Mh,z) where the luminosity ranges from the faintest 
observable magnitude (which we set to a typical ground-based sur- 
vey depth of r < 25.5), to an abs olute brightest galaxy magnitude 
Mr > -23. We use the best fit lCooravll2006 t) CLF parameters 
or the chosen fiducial values in the cases where the parameters are 
poorly constrained. This results in a realistic redshift distribution of 
the foreground galaxies with roughly 15 galaxies per square arcmin 
withz < 1.0. 

We use three simple galaxy models which we term, ellipticals, 
spirals and mix. Ellipticals are defined to have the same ellipticity 
as their parent halos. Spirals are modelled as a thick disk randomly 
misaligned with the angular momentum vector is its parent halo, 
with a mean misalignment angle o f 20° as suggested by the re- 
sults of Ivan den Bosch et alj y003) and discussed in more detail 



in H04. Note that for completeness, we also present the results, 
in section |3| for a spiral model where the disk is aligned with the 
angular momentum vector of its parent halo, as in HRH. The re- 
sults in the shear-ellipticity correlation analysis that follow in sec- 
tion|4| however, are unaffected by whether the disk is misaligned or 
not, and we therefore only present the shear-ellipticity correlation 
results from our misaligned spiral model. Our mix model popu- 
lates dark matter halos with the same elliptical or misaligned spiral 
galaxies according to the early-type and late-type CLFs of Coor^ 
( 2006). These CLFs are constrained such that they reproduce the lu- 
minosity functions and clustering properties of red and blue galax- 
ies. Our simple galaxy ellipticity models will clearly break down 
in the case of parent halos which host more than one galaxy. We 
therefore limit our study to singly occupied halos which satellite 
galaxies are unlikely to inhabit. For our mix model, this criterion 
results with roughly one third of the halos populated with elliptical 
galaxies at all redshifts. 

In this paper we determine the ellipticity parameters e-i for 
each galaxy, where Si for a perfect ellipse with axial ratio f3 at po- 
sition angle 9, is given by, 



ei 

£2 



I+/32 



cos 261 
sin 26) 



(3) 



In order to study the effect of these intrinsic galaxy shapes on weak 
lensing measurements we then define and use the alternative ellip- 
ticity parameters e^ given by 

Ei 



2R ' 



(4) 



where the responsivity Risa correction to account fo r the response 
of an ensemble of galaxies to weak lensing shear JRhodes et alJ 
2000; Bernstein & Jarvi3l2002h . Without this responsivity correc- 
tion to our ellipticity measurements, equation Q would not apply. 
For the following unweighted analysis, i? = 1 — (e^ ) — 0.89 for 
all galaxy models (although note that for simplicity in HRH and 
H04 the responsivity R was set equal to 1). The resulting a ellip- 
ticity distributions have a mean (e^) — and a width a^ = 0.33 
which is very similar to what is measured from space-based data 
(Hevmans et al.ii2006i) . 



3 INTRINSIC GALAXY ALIGNMENTS 

In this section we update the analysis of HRH and H04 by re- 
measuring the intrinsic alignment correlation (6565) from simu- 
lations. The N-body simulation used in this analysis is at a signifi- 
cantly higher resolution than the VIRGO simulations 1 Jenkins et al] 
1 1 9981) used by HRH allowing for more reliable measurements of 
halo angular momentum. It also covers a larger cosmological vol- 
ume and thus yields higher signal-to-noise measurements. The evo- 
lution of the large-scale structure is now preserved, permitting a 
direct measurement of the E and B modes of the angular intrinsic 
alignment correlation. As a result we do not have to rely on pro- 
jecting a measured real space correlation function into an angular 
correlation function using weakly constrained redshift dependent 
galaxy clustering models. The analyses of the simulations differ 
slightly in the fact tha t we now populate halos according to the 
CLF of lCooravl yOOg) instead of introducing low and high mass 
cuts to remove halos less likely to harbour a single galaxy. Our 
galaxy models are however identical, with the new addition of the 
mix model. The intrinsic ellipticity correlation rjij-ab), is given by 



■n{rab) = (et(ra)et(rb)) + (er(ra)er(rb)) 



(5) 
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Figure 1. The intrinsic galaxy ellipticity correlation riirab) as a function of 
comoving galaxy separation r^h as measured from the simulation assum- 
ing three different galaxy models; ellipticals, spirals (placed perpendicular 
to the angular momentum vector as in HRH, lower left, and misaligned, as 
in H04 and the rest of this paper, lower right) and mix. The shaded region 
shows the 95% confidence region c onstrained from measurem ents of in- 
trinsic galaxy alignments in the SDSS Mandelbaum et alJj200fil) . where the 
best-fit is over-plotted (solid). The best fit model from the misaligned spiral 
galaxies measured in H04 is shown dotted. 



where the comoving galaxy separation Tab ~ |ra — rb|, and (etet) 
and {ere,.) are the tangential and rotated correlation functions (see 
for example H04). The measured rj{rab) is shown in Figure0for 
each galaxy model, where for completeness we have also included 
the HRH spiral model where the disk is placed perpendicular to 
the angular momentum vector of the halo. The shaded region in 
each plot is the Icr confidence region for rj{rab) as constrained 
from measurements of in trinsic galaxy alignments in the SDSS 
JMandelbaum et all2006l) . The strength of the signal for the ellip- 
tical model at separations r > 3/i~^Mpc is inconsistent with the 
SDSS results at the 6a level and we can thus reject this model at 
99.9% confidence. The spiral and mix models are consistent with 
the SDSS constraints. All results are consistent with the previous 
results from HRH and H04, and now supersede those constraints. 

Following H04, we provide parameter fits to these results us- 
ing the following model 



(e,(ra)e,(rb)) 



A 



[Vab/BY 



(6) 



where the amplitude A is chosen to be a free parameter. For the el- 
liptical model, B is also allowed to be a free parameter. For all other 
models B constrained io B — l/i~^Mpc. Table Q compares the 
best fitting parameters for the tangential {etet) and rotated (erCr) 
and rotated {ere,-) ellipticity correlation functions, as well as the 
total ellipticity correlation rj{rab)- For the elliptical model, both A 
and B are free parameters. For all other models B is constrained 
to be S = l/i~^Mpc. The best fits to the results presented in 
this paper can be comp ared to the result from H04 (referred to as 
HRH* in that paper) and lMandelbaum et aljf200d) (SDSS), where 



Model 


Parameter 


{etet) 


(erCr) 


Vi^ab) 


Ellipticals 


A/10--' 


1 47+0.05 
^■*' -0.02 


r, OO+0.18 

U.d»_o.08 


1 7c: + 0.08 
^^•'3-0.02 




B 


2.94t°:°I 


n Qfi+O-io 

'j-yo_o.3o 


9 77+0.09 
^•' '-0.07 


Mix 


A/10-3 


o-45lS];J°5 


n n7+-''* 


0A7tlf, 


Spirals 


A/10-3 


n sn+005 


-0.02tH^ 


0-28l»;S° 


Misaligned 


A/10-3 


n it;+008 
U-J^o_o.02 


0.00±0;g^ 


+0.10 

^••J'-o.os 


HRH* 


A/10-3 


- 


- 


0.35 ±0.16 



SDSS 



A/10- 



0.57 ±0.72 



Table 1. The intrinsic alignment amplitude A and scale length B as defined 
in eguationlSlfor the tangential {etet), and rotated (erfir) ellipticity corre- 
lation functions, as well as the total ellipticity correlation »?(r'ai,). For the 
elliptical model, both A and B are free parameters. For all other models 
B is constrained to be B = 1 h~^Mpc. The best fits to the results pre- 
sented in this paper can b e compared to the result from H04 (HRH*) and 
iMandelbaum et all i200dh (SDSS), where we have included a scaling factor 
of 4:R? to account for the shear responsivity (equation |4} that we include 
in this analysis. The quoted errors are all Icr, except in the case of SDSS 
which corresponds to a 3(t constraint. 



we have included a scaling factor of 47?^ to account for the shear 
responsivity (equation |4} that we include in this analysis. It is in- 
teresting to note that the rotated component of the ellipticity cor- 
relation function (erCr), is found to be weaker than the tangential 
component {etet), which is in contrast to previous assumptions that 
{etet) ~ (e,.er> . 

Taking all halos between a redshift of < z < 1 we mea- 
sure the angular galaxy elli pticity correlation (es(^ a)es(6^b)), di- 
rectly from the simulations. ICrittenden et al.l i2002h show that the 
correlation functions can be decomposed into the following E- and 
B-type correlators. 



^""{0) 



^+{O)+Ci0) 



em^ 



e+w-rw 



where, in this case. 



^±(^) = (etet>e ± {erer)g, 



and 



f«-(») 



126'^ 



(7) 



(8) 



^^-W- (9) 



Following the method detailed in lPen et aU i2002h . we have 
decomposed the ellipticity correlation functions into their E and 
B modes ^^^^, assuming that {e°-e^)e = for 6 > 20 arcmin. 
Note that we treat each realization of the N-body simulations as an 
independent pointing and plot the mean and co-variance between 
the results from each realisation. The results are shown in Fig- 
ure|2|where the measured E-mode can be compared to the expected 
ACDM shear correlation from a weak lensing survey with median 
redshift Zm ~ 0.7. This redshift corresponds to the median red- 
shift of the galaxies, or the populated halos, in the simulations. For 
both the spiral and mix models, the signal is consistent with zero 
on all scales, although see H04 for a noise-free E/B mode analy- 
sis of the expected angular correlation signal from the best-fitting 
model to r]{rab), which is directly applicable to the spiral and mix 
models analysed in this work. The elliptical model exhibits strong 
E and B modes out \o 9 < 10 arcmin, although the small angle 
B-mode are consistent with zero. Figure|3]shows the E/B mode de- 
composition of the elliptical model signal when we have repeated 
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the analysis removing galaxy pairs that have a physical separation 
of rab < 10/i~^Mpc, the resulting E and B mode signal is consis- 
tent with zero on all scales, as expected. 



3.1 Intrinsic Galaxy Alignment with Voids 

We follow the analysis of lTruiillo et alj yOOg) to determine the de- 
gree of alignment between the rotation axis of galaxies located on 
the shells of cosmic voids, and the surface of those voids. We use 
a section of the simulation at fixed output time z ~ 0.325 and 
take all halos with mass Mh > 8.5 x IO^^/i^^Mq yielding a sim- 
ila r number density of the 2dFGRS and SDSS surveys analysed 
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Trujilloet alJ jgOOd) . Voids are selected using the same method 
Patiri et alj yOOg) and for each galaxy within r-void < r < 
rvoid + 4/i~^Mpc of the void centre (which has radius rvoid) we 
calculate cos{9) — r.L/(|r||L|) where L is the angular momen- 
tum of the parent halo. Figure |4] shows the resulting distribution of 
galaxies residing on the surface of voids with ^ 9 ^ 90 degrees. 
The errors come from a bootstrap analysis, and the solid line is the 
average result when the halo angular momenta are randomised. We 
find a distribution that is consistent with the null hypothesis that the 
angular momentum vector of halos near voids is randomly oriented 
with respect to the centre of voids. For our spiral galaxy model, 
where galaxy orientation is determined by the angular momenta of 
the parent dark matter halo, this result implies that no significant 
correlation exists between the orientation of galaxies and the cen- 
tre of nearby voids in our simulations. This is in contrast to the 
[Truiillo et al. (2006) detection results, which are over-plotted on 
figure|4|that can reject this null hypothesis at the 99.7% confid ence 
level. These results are well fitted by the analytic prediction of lLed 
i2004) . (shown dot-dashed). We discuss this difference further in 
section|5| 



4 SHEAR-ELLIPTICITY CORRELATIONS 

In this section we present shear-ellipticity correlation functions and 
their E/B mode decomposition for the three different galaxy mod- 
els, as described in section ITTI Figure |5| shows the measured cor- 
relation {'ye) between foreground galaxy ellipticity e (where the 
foreground galaxies have redshifts < z < 1) and the gravita- 
tional shear 7 experienced by a set of background source galaxies 
at 2. ~ 1 where 



{je)e = (7t(6'a)et(6»b)> + (7r(6'a)e,(eb)> 



(10) 



Note that we treat each realization of the N-body simulations as 
an independent pointing and plot the mean and variance between 
the results from each realisation. This result represents the level of 
contamination that one would expect from each galaxy model for a 
deep lensing survey with a median redshif t Zm ^ 1- These results 
show an anti-co rrelation, as predicted bvlHirata& Seliald i2004) 
and measured bv lMandelbaum et alJ i200(^. The strongest signal is 
found from the elliptical model, where all occupied halos contain 
a galaxy whose shape is the same as its parent dark matter halo. 
The spiral model shows a correlation that is consistent with zero 
on all scales, and the mix model shows a weak but significant anti- 
correlation on scales 9 < 20 arcmin. For reference, the solid theory 
curve shows the the shear-shear correlation function (multiplied by 
— 1) for a deep Zs = 1 survey using the adopted ACDM cosmology. 
Note that the shear-ellipticity correlation for the mix model has a 
similar magnitude to the models tested in King ( 2005|). In contrast 
however, we find signal that has power out to angular scales that 
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Figure 2. The E mode (circle) and B mode (cross) decomposition of the 
angular intrinsic galaxy ellipticity congelation (eje^)^ as a function of an- 
gular galaxy separation on the sky, 6. This is measured from the simulation 
assuming three different galaxy models; ellipticals (upper panel), spirals 
(middle panel) and mix (lower panel). The resulting E-mode can be com- 
pared to the expected sheai' correlation (solid line) from a weak lensing 
survey with median redshift Zm = 0.7, which corresponds to the median 
redshift of the galaxies in the simulations. For clarity the B-mode points are 
slightly offset to the right. 
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Figure 3. The measured angular intrinsic galaxy ellipticity congelation 
(eje^)^ as a function of angular galaxy separation on the sky, 8 when 
galaxy pairs that have a physical separation of Tab < 10/i~^Mpc are re- 
moved from the analysis. This is measured from the simulation assuming 
the elliptical galaxy model (right) and mix galaxy model (left). The shear 
correlation from a weak lensing survey with median redshift Zm = 0.7, 
which coiTesponds to the median redshift of the galaxies in the simulations, 
is shown for reference (solid line). 
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Figure 4. The distribution of galaxies residing near to a void as a function of 
their degree of alignment between the rotation axis of the galaxies and the 
surface of the void (circles). For a perfectly aUgned set of galaxies we would 
expect to see a peak at 6 = 90° . Instead our result is consistent with the 
null hypothesis that galaxies ai'e randomly oriented with all values of cos(9) 
equally likely. The errors come from a bootstrap analysis, and the solid 
line is the average result when the halo angular momenta ar e randomised. 
Our re sults can be compared with the observational results of lTruiillo et alj 
12009) (crosses) and the best fitting analytical prediction to these results 
from lLed i2004) . shown dot-dashed. 



are sig nificantly larger than the 1 arcmin cut-off considered in lKina 
i2005 1. 

Table |2| details the level of contamination to a weak lensing 
measurement by a shear-ellipticity correlation from each galaxy 
model, for both a deep Zs ^ 1 survey, whose results are plotted 
in figure |5| and medium deep Zg ~ 0.5 survey. The mix model 
shows a level of shear-ellipticity correlation which is —6 ± 1% of 
the lensing signal at an angular scale 6 — 1 arcmin for the deep 
Zs ~ 1 survey. 

Figure |6| shows the shear-ellipticity correlation function {'je) 
for the mix galaxy model for two samples of galaxies split by 
mass. The strongest shear-ellipticity correlation is seen for the most 
massive galaxy sample where Mh > 1.0 x lO^^/i^^M©. This 
is expected, as the strength of the lensing shear is proportional 
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Figure 5. The shear-ellipticity correlation function (76)9 for the elliptical 
(circles), spiral (triangle) and mix (filled) galaxy models. In this result all 
source galaxies (that experience a lensing shear 7) are placed at a source 
redshift Zs = 1. The foreground galaxies, with intrinsic ellipticity e, pop- 
ulate the redshift space < z < 1. These results are tabulated in table I2I 
along with the results for a survey were all source galaxies are placed at 
Zs = 0.5. For reference, the solid theory curve shows the shear-shear cor- 
relation function (multiplied by —1) for a deep 2s = 1 survey using the 
adopted ACDM cosmology. 
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-67 ± 5 % 

8±5% 
-24 ± 4 % 


-43 ± 5 % 

3±5% 
-14±5% 


-82 ± 8 % 
-10±8% 
-31 ± 8 % 



Table 2. The percentage of contamination to a measurement of the shear- 
sheai' correlation function by a shear-ellipticity coiTelation from each galaxy 
model 3X9 = 0.5, 1.0, 10.0 arcmin. Results are given for both a deep Zs ~ 
1 lensing survey (upper rows), whose results are plotted in figure |5] and 
medium deep lensing Zs ~ 0.5 survey (lower rows). 



to the mass of the lens. I t is also reminiscent of what in seen in 
iMandelbaum et al.l i2006l) . where the strongest signal is measured 
from the most luminous and hence the most massive galaxy sam- 
ple. This result allows us to briefly discuss the impact of multiply- 
occupied halos, which we have removed from our analysis as a 
result of our uncertainty in how to connect the properties of the par- 
ent dark matter halo to the luminous galaxy members. One could 
consider the case where the central galaxy is represented by our el- 
liptical model and the satellite galaxies are randomly oriented. The 
multiply occupied halos are typically the more massive halos, and 
we would therefore expect to see a strong signal from the central 
galaxies, as seen with the most massive galaxy sample shown in fig- 
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Figure 6. The shear-ellipticity correlation function {-fe)g for the mix 
galaxy model for two samples of galaxies split by mass; where M^ > 
1.0 X IQI^/i-^Mq (open) and where M^ < 1.0 X IO^^/i-IMq (filled). 
Note that all source galaxies (that experience a lensing shear 7) are placed 
at a source redshift 2,, = 1 and that for reference, the solid theory curve 
shows the the shear-shear coiTelation function (multiplied by — 1) for a deep 
Zs = 1 survey using the adopted ACDM cosmology. 



ure|6| Randomly oriented satellite galaxies would not contribute to 
the shear-ellipticity correlation signal and would thus dilute effect 
of the central galaxy. The reality of this scenario is however diffi- 
cult to quantify, and indeed there is evidence that suggests satellite 
galaxies are not randomly orient ed with a tendency to align with th e 
major axis of the central galaxy terainerd''2005VYane et al."200(?). 
We therefore limit our discussion, in the rest of this paper to singly- 
occupied halos, deferring a quantitative analysis that includes satel- 
lite galaxies for a future paper analysing higher resolution N-body 
simulations. .^^^^^^^^^ ^_^ 

^^_ Following Crittend en et alJ J2003) and the method detailed in 
|Pen et al.l200Z) . we decompose the measured shear-ellipticity cor- 
relation functions into their E and B modes, equation Q assuming 
that {'ye)g — for 9 > 20 arcmin. Figure shows the result- 
ing E and B mode signals for all three galaxy models. As expected 
the E and B modes from the spiral model are consistent with zero. 
For the mix and elliptical models, the B-modes are negative but 
weaker than the corresponding E-modes. The B-modes are how- 
ever fairly constant on scales 9 < 15 arcmin and would therefore 
be detectable if all B-modes from other non-lensing sources, such 
as residual point-spread function distortions, are negligible. 



4.1 Redshift dependence 

In this section we measure the dependence of the shear-ellipticity 
correlation signal as a function of the foreground galaxy redshift. 
We present the results for the elliptical model as the strength of the 
signal allows us effectively to probe the redshift dependence in the 
presence of increased noise from the redshift binning. The trends 
seen with the elliptical model are also seen with the mix model. 
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Figure 7. The E mode (circle) and B mode (cross) of the shear-ellipticity 
correlation function for the elliptical (upper panel), spiral (middle panel) 
and mix (lower panel) galaxy models. For clarity the B mode points are 
slightly offset to the right. 



although at lower significance, and the shear-ellipticity correlation 
for the spiral model remains consistent with zero for all redshift 
slices. 

Figure |8| shows the redshift dependent shear-ellipticity corre- 
lation for 5 different galaxy redshift slices; 0.0 — 0.4, 0.4 — 0.65, 
0.65 - 0.75, 0.75 - 0.85, 0.85 - 0.95. This result demonstrates 
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Figure 9. 2ct constraints on the amplitude A and scale dependence do 
(equation I12i of the shear-ellipticity coirelation assuming the elliptical 
galaxy model. Parameter constraints are determined from the 5 different 
galaxy redshift slices shown in figure|8]where 7 is measured at Zs = 1.0, 
plus 2 additional low redshift slices where 7 is measured at Zg = 0.5. 
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Figure 8. The redshift dependence of the shear-ellipticity coirelation as 
a function of galaxy redshift slice zi. The upper panel shows the mea- 
sured signal for the following galaxy redshifts; 0.0 < 2; < 0.4 (circle), 
0.4 < zi < 0.65 (solid), 0.65 < zi < 0.75 (square), 0.75 < 2, < 0.85 
(hexagon), 0.85 < 2; < 0.95 (triangle) . The dashed line shows the ellip- 
tical measurement in Figurelslfrom the full galaxy sample 0.0 < 2j < 1.0. 
The lower panel shows the same signal with the lensing efficiency E fac- 
tored out and normalized to the lensing efficiency N of the mean galaxy 
redshift of the full galaxy sample where 2; = 0.68 



that the most significant anti-correlation is found for the low red- 
shift (z < 0.65) galaxy samples, which are the most efficient for 
lensing source galaxies at Zs = 1. The lensing efficiency E of the 
lens-source pair is defined to be 



E = 



DiDl 



where the distances D are all comoving angular diameter distances 
to the lens D; , to the source Ds , and between the lens and source 
Dia. We might expect the shear-ellipticity correlation signal to 
scale with the lensing efficiency and to test this hypothesis, the 
lower panel of figure |8| shows the shear-ellipticity correlation sig- 
nal with the lensing efficiency factored out; we divide by E for the 
median lens redshift in each redshift slice, and multiply by a nor- 
malising factor A'' set equal to the lensing efficiency E of the mean 
redshift lens for the whole halo sample where zi = 0.68. The result 
shows that the redshift dependent shear-ellipticity correlation mea- 
surements are very consistent with one another once the lensing 
efficiency is factored out. This then motivates the following param- 
eterisation of the redshift evolution of the shear-ellipticity correla- 
tion. 



{'y{z,)e{zi))e = E- 



A 



+ 00 



(12) 



(11) 



where the amplitude A and scale dependence 6*0 are free parame- 
ters. E is the lensing efficiency of the lens-source (76) pair (equa- 
tion lllk and zi is set equal to the median redshift of each galaxy 
redshift slice. Note that this parameter isation is sim ilar to the ex- 
ponential parameterisation proposed bv lKingH2005l) . but this func- 
tional form is found to produce a better fit to the results. Figure |9| 
shows the two parameter 2a constraints on the amplitude A and 
scale dependence 9o as measured independently from the 5 differ- 
ent galaxy redshift slices shown in figure|8|where 7 is measured at 
Zs — 1.0, plus 2 additional low redshift slices where 7 is measured 
at Zs — 0.5. The consistent results for each galaxy redshift slice 
suggest that the evolution of the shear-ellipticity correlation can be 
safely modelled through the lensing efficiency of the lens-source 
pair. 

Using this parameterisation (eQuation ll2t we place constraints 
on the amplitude A and scale dependence 60 of the shear-ellipticity 
correlation measured jointly from the redshift dependent results. 
Figure fTol shows the resulting two-parameter 1 and 2a constraints 
for the elliptical model (upper) and mix model (lower). The best 
fitting results for the elliptical model are A = (—5.60 ± 0.18) x 
W^hMpc"^ and ^0 ~ 1-83 ± 0.13 arcmin. For the mix model 
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Figure 10. 1 and 2cr constraints on the amplitude A and scale dependence 
60 teauation ll2l of the shear-ellipticity correlation (76) , assuming the el- 
liptical galaxy model (upper) and mix galaxy model (lower). For clarity, 
only the la constraints are shown for the tangential {"ftet) and rotated 
{"fr^r) shear-ellipticity correlations, where {'^re-r) is the weakest signal. 
The best-fit parameters are maiked with a stai'. 



A = (-1.29±0.16)xlO"^/iMpc-^ and6lo = 0.93J:[;:g arcmin. 
The Icr errors on A are marginalised over 60 and vice versa. The 
constraints measured for the tangential {'ytet) and rotated {'yrer) 
shear-ellipticity correlations are also shown in figure [Tol but for 
clarity only the la constraints are presented. Note that the (7,6^) 
shear-ellipticity correlation is always weaker than the {'jtet) shear- 
ellipticity correlation. 



5 DISCUSSION 

The results that have been p r esented in this paper and o thers (e.g 
H04. lHirata & SeliakI i2004l) : lMandelbaum et alj J200d) ). demon- 
strate that if a relation exists between the shape of galaxies and the 
properties of their parent dark matter halos then there is the poten- 
tial for contamination to the measurement of weak gravitational 
lensing by large-scale structure. The conclusion that one should 
draw from these works is that it is no longer acceptable to assume 
that all the intrinsic contributions to the observed galaxy elliptic- 
ity correlation function (i.e. terms that include Ca in equation |2j 



are negligible. We must therefore seek to derive and apply optimal 
methods to remove these sources of contamination in the analyses 
of the next generation of weak lensing surveys. In this paper we 
have sought to answer the question, what level of contamination 
might one expect, and our conclusions are clearly dependent on 
our choice of galaxy model; ellipticals, spirals or mix. 

A wealth of observational evidence allows us immediately 
to rule out the elliptical model on several grounds. The mea- 
surement of the intrinsic alignment signal ri{r) from the SDSS 
iMandelbaum et alJl2006l) is inconsistent with the elliptical model 
signal at the 99.9% level, as shown in FigureQ Similar constraints 
could also be placed from the intrinsic alig nment signal measured 
in the SuperCOSMOS survey (Brown et al. 2003). Figure|2|shows 
the expected angular ellipticity correlation in the absence of weak 
lensing for a survey with median redshift Zm = 0.7. This medium 
depth survey is similar to the Red sequence Cluster Survey (RCS) 
[Hoekstra et al. (2002) and the CTIO survey (Jarvisetal. 2003J). 
Weak lensing analyses of these two surveys constrain the amplitude 
of the matter power spectrum to be ug ~ 0.8 ± 0. 15 for Qm =0.3, 
which is consistent with th e results from the an alysis of the cos- 
mic microwave background (Soerael et al. 2006). If more than half 
of the RCS and CTIO measured lensing signal came from intrinsic 
alignments, as suggested by the results of the elliptical model, the 
'true' value of as would be then be very low and inconsistent with 
other observational constraints. The COMBO-17 survey has been 
analysed using redshift weighting to remove the intrinsic contribu- 
tion from close galaxy pairs (H04) and only a small difference was 
noted in the resulting signal, placing constraints on the intrinsic 
alignment contamination that would also reject the elliptical model 
at the 99% confidence level. Lastly, if the shear-ellipticity correla- 
tions were as strongly anti-correlated as expected from an ellipti- 
cal model for median-deep surveys (see Table|2j we would expect 
to measure no weak lensing signal in the RCS and CTIO surveys, 
which is not the case. All of the evidence therefore suggests that 
not every galaxy has the same shape as its parent halo, and we can 
therefore rule out the 100% elliptical model. 

The results from th e spiral model are consistent with the 
iMandelbaum et al] yOOg) SDSS constraints on intrinsic galaxy 
alignment. It predicts a very low intrinsic alignment correlation and 
low shear-ellipticity correlation that is consistent with zero on all 
scales. It is important however to take note o f the SDSS results 
for their 'L6' massive elliptical galaxy sample JMandelbaum et alJ 
2006). This sample of massive ellipticals exhibit a strong correla- 
tion between their galaxy shape and the surrounding density field. 
This observation has been propagated to high redshift to predict a 
measured shear-ellipticity anti-correlation that is ~ 100% ± 50% 
of the magnitude of the expected shear-shear correlation signal for 
a survey with Zm = 1- The amplitude of the 'L6' shear-ellipticity 
correlation is similar to the measured shear-ellipticity correlation 
signal from our elliptical model, but is inconsistent at the 4cr level 
with the measured shear-ellipticity correlation signal from our spi- 
ral model. In order to reproduce the strong SDSS signal from a sub- 
sample of galaxies in our simulations we therefore require a small 
population of elliptical galaxies and can thus rule out the 100% 
spiral model. 

We have modelled a mixed population of spiral and elliptical 
galaxies using the early-type/late-type halo occupation distribution 
of Coorav ( 2006) which reproduces the observed luminosity func- 
tions and clustering properties of red and blue galaxies. The results 
from this mix model are consistent with the all the observational 
results ofjMandelbaum et al. (2006) It predicts low intrinsic align- 
ment correlation, weak, though significant, shear-ellipticity corre- 
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lations, and negative B-modes. B-modes from systematics associ- 
ated with the observations (for example residual distortions from 
the optical point-spread function) would typically produce positive 
B-modes and hence one would not expect to be able to detect the 
weakly negative B-mode signature of the mix model in current data. 
The effect of the shear-ellipticity correlation is to reduce the 
measured shear-shear correlation below the amplitude of the true 
cosmological lensing signal. If uncorrected for, this would result in 
an underestimate of cts from weak lensing measurements. The de- 
gree of underestimation would be survey depth dependent as from 
Table |2| we expect a medium-depth survey to suffer ~ 20% con- 
tamination and a deep survey to suffer ^ 10% contamination. A 
comparison of the results from different depth surveys should there- 
fore, in principle, be able to constrain models of shear-ellipticity 
correlations as one would expect to measure lower values of ag, 
from shallower surveys. Whilst it is tempting to discuss the trend 
in the literature of somewhat higher ag, measurements from deeper 
lensi ng surveys, the unce rtainties in the source redshift distribu- 
tion iVan Waerbeke_et_alj [20()6). sampling variance and shear cal- 
ibration errors faevmans et alj 2006) means that we are unable to 
make any firm conclusions with current lensing data to confirm or 
rule out the mix galaxy model. 

5.1 Diagnostics for shear-elllptlclty contamination 

For the next generation of lensing surveys, the technical sys- 
tematics associated with the detection and measurement of weak 
lensing have a promising chance of being negligible. A signifi- 
cant detection of negative B-modes would then suggest contam- 
in ation from shear-ell ipticity correlations. In addition, the results 
of iMandelbaum et alJ i200a) and the comparison of the different 
models used in this analysis demonstrate that shear-ellipticity cor- 
relation is likely to be galaxy morphology dependent. For future 
multi-colour surveys with photometric redshifts one will be able 
effectively to separate red-sequence elliptical galaxies from the 
galaxy sample using a simple colour cut (see for example lBell et al] 
yOMD). With large area surveys one could then analyse the data 
with and without the elliptical galaxy population. A difference be- 
tween the two measurements would suggest contamination from 
shear-ellipticity correlations. If that is the case, a zero B-mode mea- 
surement from an analysis where all elliptical-type galaxies have 
been rejected might well serve as a good indication that the con- 
tamination has been removed. 

In section l4m we investigated the redshift dependence of the 
shear-ellipticity correlation (7c, e^) where Za < Z0. The results, 
shown in figure |8| demonstrated that the evolution of the signal is 
proportional to the lensing efficiency E (eguation llH of the lens- 
source pair, where, for the shear-ellipticity correlation {'-faep) the 
source is at Za and the lens is at zp. The model that provides the 
most realistic description appears to be the mix model for which we 
have constrained the amplitude of the {-ye) signal to ~ 10% accu- 
racy. We have calculated the level of contamination arising from a 
mixed galaxy population for a deep Zm ~ 1 survey analysed using 
the shear correlation function (77), and found it to be at the level 
of —6 ± 1% and — 10 ± 2% at angular scales I and 10 arcmin re- 
spectively. Using the (76) model (eauation ll2> with the best fitting 
parameters to our results, one could now potentially remove this 
contaminating signal to an accuracy of better than 1 % of the lens- 
ing signal. We would urge caution with applying such a method 
however, as our results in section ITTI clearlv indicate the limita- 
tions of the simple models that we have used to link the properties 
of galaxies with their parent dark matter halos. We would there- 



fore recommend using the (76) model (equation 1 12t that we have 
provided as a guide to the level of contamination, which should be 
self-consistently determined from the lensing data using the result 
th at the contaminat ion scales with lensing efficiency. As proposed 
bv lHirata & SeliakI J20 04'). measure ments of cross-correlati on to- 
mography {7a 7/3) (see for example iTakada & Whitd i2004h ') will 
be a useful diagnostic tool, as the true cosmological lensing sig- 
nal scales differently with redshift with the lensing efficiency in- 
tegrated over all foreground lenses. One would therefore able to 
project out the diff erent redshif t dependences of the two signals 
as demonstrated bv lKina y003) and simultaneously constrain the 
cosmological lensing signal and the contaminating shear ellipticity 
correlations. 



6 CONCLUSION 

We have investigated potential sources of contamination to the 
measurement of weak lensing by large-scale structure from the in- 
trinsic alignment of galaxies, in section|3| and from the presence of 
shear-ellipticity correlations between foreground galaxy shapes and 
the weak lensing shear experienced by background galaxies, in sec- 
tion|4] These two sources of physical contamination behave rather 
differently. The intrinsic alignment contamination arises from close 
galaxy pairs, and can thus be eliminated by downweighting those 
galax y pairs which are at clos e separations JHevmans & Heavens! 
2003l lKing & Schneidedl2()0a) . or by measuring cross-correlation 
tomography between different redshift bins. Intrinsic ellipticity- 
ellipticity correlations are thus not a problem if photometric red- 
shifts are available. In contrast, the shear-ellipticity correlation 
arises from galaxy pairs which are at large physical separations, 
with the strongest correlations arising for lens-source pairs where 
the lensing efficiency E (eguation llH is at its strongest. This could 
potentially be directly eliminated through morphological galaxy se- 
lection, or from measuring weak lensing shear correlation within 
thin redshift slices (Kins 2005). All possible techniques to remove 
these sources of contamination will require good photometric red- 
shift measurements and acquiring deep multi-colour data should be 
made a top priority for all future lensing surveys. 

As discussed in section|5|the conclusions are dependent on the 
choice of model that we use to connect galaxy shape to the proper- 
ties of its parent dark matter halo. General trends are that for deeper 
surveys, where the cosmological lensing signal is at its strongest, 
the degree of contamination is decreased. Eliminating the elliptical 
model through comparison to several observation results allows us 
conclude that for medium deep surveys Zm > 0.5, intrinsic align- 
ment contamination contributes only to the small angular scales 
6 < 1.5 arcmin, and we refer the reader to H04 for a more de- 
tailed discussion on the implications of intrinsic galaxy alignments 
for future lensing surveys. The shear-ellipticity correlation is con- 
sistent with zero for our spiral model. To reproduce the results of 
IMandelbaum et alj J 20 0o) however it is necessary to include a pop- 
ulation of elliptical galaxies which introduce contamination to the 
shear lensing signal at the ~ 10% level for a deep Zm ~ 1 survey 
and ~ 20% level for a medium deep Zm ~ 0.5 survey (see table|2j. 

The galaxy models that we have applied to the simulations 
are indeed clear simplifications of reality as they consider only the 
properties of the parent halo, not the halo environment or the tidal 
density field directly. The non-detection of any alignment between 
a galaxy's spin axis and the surface of voids in our simulations, 
in section Ism i s in some tension w ith the analysis of the 2dF and 
SDSS surveys JTruiillo et alj|200a) . If the observation is robust, it 
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is therefore likely that it is not only the properties of the parent 
halo that impact on galaxy shape, but also the surrounding envi- 
ronment to some extent. In the rare cases of galaxies near voids or 
clusters, we might well expect the environment to have a stronger 
effect on galaxy shape than its parent dark matter halo and it will 
be important to investigate this further with full hydro-dynam ical 
simulations on cosmological scales (see for example Navarro e t alj 
|2004 . 

The important point to take away from this analysis is that our 
numerical study supports the notion that the contamination from 
shear-ellipticity correlation scales with lensing efficiency. This sim- 
ple scaling with redshift offers possibilities of determining the cos- 
mic shear signal and the contamination simultaneously, either by 
using the models constrained in section l4m for the behaviour of the 
contamination ("cf lKinal2005h . or by making fewer assumptions and 
parameterising the contamination. Because of the redshift scaling, 
only a single function needs to be parameterised. Cosmological pa- 
rameter estimation would then proceed by marginalising over the 
nuisance parameters describing the contamination. Within the con- 
text of the model considered here, the (76) contamination can be 
estimated to an accuracy of 3% and 10% for the elliptical and mix 
models respectively. 

A more detailed study of the redshift dependence of shear- 
ellipticity correlation will, however, be required for the more am- 
bitious future lensing experiments which aim to determine w{z). 
Such studies will need to include the satellite galaxy population 
that is missing from this analysis in addition to an improved mixed 
galaxy population model that takes into account an observation- 
ally constrained redshift evolution in the ratio of elliptical to spiral 
galaxies for different halo masses. This will require the develop- 
ment of better models for galaxy formation in dark matter simu- 
lations to improve our understanding of the connection between 
galaxies and their dark matter halos and the resulting impact on 
weak lensing studies. In the meanwhile we should rely on prob- 
ing this phenomenon with existing data sets, in order to develop 
suitable diagnostics and analysis tools that will not subject weak 
lensing measurements to the kinds of contamination investigated in 
this paper thus allowing weak lensing to reach its full potential as a 
high precision dark matter and dark energy probe. 
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